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In previous papersi'i we showed that a division of the /-sheU into two subsystems gives a better 
understanding of cohesive properties as weU the general behavior of lanthanide systems. In this 
paper we present numerical computations using the suggested method. We show that the picture 
is consistent with most experimental data, e.g. the equilibrium volume and electronic structure in 
general. Compared to standard energy band calculations, or to calculations based on the self inter- 
action correction and LDA+U, in the present approach /-(non-/)-mixing interaction is decreased 
by spectral weights of the many-body states of the /-ion are involved. 



I. INTRODUCTION 

Strongly correlated electrons show a rich variety of physical properties, e.g. heavy Fermion behavior, partial local- 
ization and unconventional superconductivity. Normal ab initio electronic structure calculations performed for these 
systems have proven to be unable to capture the essential features of these systems. The reason for this is mainly 
that the LDA underestimates the strong correlations caused by the local part of the Coulomb interaction among the 
electrons. The standard model for the lanthanides corrects for this via imposing a condition on the /-electrons which 
is not present in the standard LDA DFT scheme: i) /-electrons should be treated as atomic like and localized; ii) 
the /-shell contains an integral number, n, of electrons. However, it has been argueda, that already the cohesive 
properties signal a dev,iation from the standard model. 

In a previous papera, we presented an idea how normal band structure calculations could be extended to treat 
systems with strongly correlated electrons, that is an improvement on normal methods such as LDA+U, SIC or 
treating the f-electrons as atomic like. Indeed it was demonstr|ated that the cohesive properties, such as e.g. equilibrium 
volume is better described with this method. The calculationo has been performed within the linear muffin-tin orbital 
method within the spherical approximation (LMTO-ASA). The comparison with the results of normal LMTO-ASA 
calculation has shown that the latter contains two type of errors: one comes from the spherical approximation and 
another from the underestimated electron correlations. Here we have implemented these many-body corrections in 
a full-potential all-electron relativistic energy band method, in order to obtain a description as accurate as possible, 
and we have tested our theory for the light lanthanides. These elements are an excellent testing ground for any theory 
that tries to incorporate electron correlations, since so much is known experimentally about them. For this reason, we 
have chosen to test our theory on these materials, before addressing more involved systems such as the heavy Fermion 
materials, UPta and CeCu2Si2 and possibly the 5-phase of Pu. 

The paper is organized as follows. In Section II we derive the equations that form the ground for the calculations, 
in Section III we describe the details of the numerical implementation, and Section IV contains the results. Finally 
we conclude in Section V. 

II. DERIVATION OF THE KOHN-SHAM LIKE EQUATIONS 

The main motivations of this work are to derive expressions for corrections from strong electron correlations to 
standard methods of ab initio methods. We also want to make an attempt to understand, from a many-body 
theory point of view, the mechanism of localization of f-electrons, that may be described phenomenologically by 
self- interaction correction (SIC), and the so called LDA-I-U method. As we shall see the method allows us to capture 
features of correlated systems such as cgpthe correlation narrowing of bands. The full description of our approach, 
as well as notations, will be given in ref.l3, and we present here only some of the theoretical details, how we have 
implemented the theory numerically as well as the results of some of aur calculations. 

We formulate the problem following the ideology presented in reffl, namely, that the one-to-one correspondence 
between the ARF (the approximation of renormalized Fermionso) series and the one from weak-coupling perturbation 
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theory allows to use the same analytical form for the exchange-correlation potential. Moreover, our theory includes 
renormalized mixing, hopping, overlap matrixes and charge density, by the spectral weights of the many-electron ion 
states that we involve. The ARF can be viewed as the Hubbard-I approximation which is extended to the case of 
full Coulomb interaction. This approximation is the simplest possible, but since it is exact in the atomic limit it is 
expected to give reasonable results in the case of very strong correlations. 

Now let us establish a connection between the secular equations solved in electronic structure calculations (ESC) 
and Green function methods for strongly correlated electrons. In order to do this we need to define the Hamiltonians 
for the two problems. First, when solving the Kohn-Sham equation we have to solve a secular equation, which looks 
like, 

{<^OjL.nL, - h^jL.nLjXnLi = 0, (1) 

where OjLjiLi is the overlap matrix between the cite-centered basis-functions, XjiLi', j is the site index and L is a 
composite index L — [I, mi, a). ejL is the energy for the orbital, and /i^^ is the Hamiltonian matrix element which 
will be defined below. Let us first consider the /-electrons as localized. We, therefore write the Hamiltonian for the 
band structure method with the /-electrons separated from the other conduction electrons. Denoting the operators 
for /- and conduction-electrons with / and c, respectively, the ESC Hamiltonian is 

k,Li,L2 jii^i 

+ "V t"^"^ ft f ■ 

3132 J 3l'^l-l 321^2 

+ E ^3lLlj2l^2[^^Likfj2iy2~^fj2l^2^Lik]' (2) 
k.Li J2l'2 

where the matrix elements comes from the expectation value of the Hamiltonian with respect to the basis functions 

/ii„L.(k) = (XL,(k)|H|xL.(k)). (3) 

The index L is used solely for conduction electrons. We denote the /-/ hopping by t and the mixing between c and 
/ by V. Combined indices for the /'s are denoted with v. The operators fulfill the commutation relation 

where a is either a c— or an /-operator. The Hamiltonian, Eq.(^), is a DPT Hamiltonian using either LDA or any 
other approximation for the exchange-correlation potential, i.e. 

where ^''^^(^(r)) is the exchange correlation potential and Ven is the interaction of electrons with nuclei. From these 
definitions we can define Green functions which we can use to make a comparison to the case for strongly correlated 
electrons. 

Obviously, the same problem can be written in terms of Green functions (GFs). Denoting the Fermionic Green 
functions F and separating the c— and /-electrons we can write the GF, F, in super matrix form, 

/ F-;,(k,z^) Fl{likM\ 

* K^'^ ^3lLi^32L2) ^'^ ^ 3ilJ 32L2) ) 

where T is the chronological ordering operator. Further, we denote the time Fourier transform by a small to and k- 
space Fourier transform with a small k. Taking into account the commutation relations, Eq.(^ in this representation, 
we obtain the equation of motion for the Green functions (the k in 0(k) is omitted for convenience): 

[5ll21^ - Ol\hl^^^ - Ollv:^^^] -\Olly.2.. + OllVi^,,,\ \ ( ^^£f^,(k, za;) F'^^^^V, tu;) \ 



oi\, Oil, 
o-l o-l, 



(7) 
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Thus, if we multiply both sides of Eq.(0) by the overlap matrix, then, the inverse of the GF coincides with the bracket 
in Eq. (|l|) , and the corresponding secular determinant gives a band structure coinciding with Eq. . Let us now move 
to the problem of strongly correlated electrons. We will work with the multiple-orbital periodical Anderson model 
(MO-PAM) in terms of Hubbard operators. The latter are defined as projection (transition) operators acting on the 
many-body states of an ion, i, in the position Ri, as. 



(8) 



Here |r„) is an n-particle state constructed out of orbitals from the ESC, |r„) = /1/2 ■ • ■ /n|0). To simplify the 
notations we denote the Fermi-like transitions [r,r'] with Roman letters a,5, . . .. The commutator relation for the 
X-operators are written in terms of structure constants e of the algebra 



(9) 



where Z^, denotes a Bose-like transition, changing the particle number with an even number. The Fermion /-operators 
can be expressed in terms of AT-operators (or vice versa), i.e.. 



(10) 



r,r' 



The model for the strongly correlated electrons, the multiple-orbital periodical Anderson model (MO-PAM) (Eq. 
with hopping t^O), is expressed in terms of A-operators using the relation in Eq. (p^) , 



k,Li ja 



-f- c.c 



(11) 



k,Li jV2,a 



where the band-structure for the c electrons is ekLi and = £°(/j,y2)"(/,v2)"- Within the ARE we are dealing with 
the same effective Hamiltonian as Eq.(H) with the exception that instead of the term ~ have to write the 

term describing the many-electron states of an ion. 



(12) 



and /-operators should be rewritten in terms of Hubbard operators (see relO). Then the effective Hamiltonian acquires 
the form of the multi-orbital periodical Anderson model (MO-PAM). The equations for the Green functions for the 
multiple-orbital periodical Anderson model (MO-PAM) in the Hubbard-I approximation within a non-orthogonal 
basis set have the following form. 



'^jL,jiLi"-jiLi,]2L2 



J2L2 J 



SjjjLL2idt 

= z5(<-0({c,x,4}) 

+c^jz^,,j^j''^.,,2L2(i^)'^4'(4)G;l,,'(M'), 



(13) 



[zdt-^a]GfZj'=^&{t-mX'^,ri}}) 

I (H-i /,c + n^^ V* 

^ [ jti,jlLl"'jlLi,j2L2 '^jfJ.l,j2L2 

<fltef{Z^)GT2Lrm 

+OjLnLV,L,,^.MlM''izl){U^^^^^ (14) 
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where ij denotes either a c or a X-operator, is the energy of the transition a (e.g. if a = r„_i,r„ then Aa = 
Er„ — £^r„_i). We simphfy the model further neglecting the energy differences between different transitions in the 
lower Hubbard sub-band, i.e., between n- and {n — l)-electron states. Then all these transitions have the same energy, 
Ai. The same will be assumed about the upper transitions, between {n + 1)- and n-electron states; they then have 
the energy A2 . We assume that the lower transition lies so much lower in energy than the bottom of all conduction 
bands that the overlap integrals as well as mixing interaction of conduction electrons with these transitions are equal 
to zero. Therefore, for a ~ [r„_i^r„], the GFs Gfa''j' — and 

^]a,j'a' = iuj - Aa' ^^^^ 

This GF gives the equation for the population numbers 

iVr„ = (^r„_i + A^rJ/(Aa), (16) 

where /(Aq.) is Fermi function for the transition A^. Since /3Aq, ^ 1, A^r„_i — exp (— /3Aa)iVr„ . The lower population 
numbers are negligible and only Nr^ and -/Vr„+i have nonzero values, while Eq.(|l^) becomes an identity. This allows 
us to formulate the system of equations for GFs which contains only upper transitions. After Fourier transformations 
t — t' ^ iuj and Rj — Rj' k, the system of equations acquires the following form: 



+ [OZi(k)(/p)"i?ref + OZi,(k)X^L,^(k)(^)^]Gf'''(k,^c.) 
+OZ^^(k)(^J'^y;i,(k)(/t)'^ef(Z«)G^f(k,zc.), (17) 

[loj - A.lGf' (k, iLo) = i5{t - t')({x;, ,7],})k 

+0~l^S'^)VL,^{mlf^f{Z'){f,YG^^\\.M- (18) 

Let us now try to identify if we, from the general equations, can establish some connection to ESC. In ESC calculations 
n orbitals are localized and (for /-systems) 14— n orbitals are delocalized. This picture naturally arises in our approach 
within the orbital representation, where the wave function is taken in the simple form |r„) = /^j/^^.../^^ |0) = |0). 
Since in this case the occupied orbitals are fixed, there are only one such state |^ Then, the states with [n + 1) 
electrons arise when one of the orbitals ip^^, ipv2, Vvu-^ j which are empty in the state |0), is occupied in the state 
|r„_|_i) = /jj|0). Any of the states |r„+i) contains only one /-orbital of the type u (which is delocalized due to its 
overlap and mixing with the conduction-electron orbitals), each particular transition a : |r„+i) — > |r„), or, more 
explicitly, fl\Q) |0), described by the Hubbard operator X'^, in this representation can be relabeled by X'' since it 
fills the orbital v. Thus, one can indeed establish a one-to-one correspondence between filling unoccupied /-orbitals 
and the many-electron transitions. But this is possible only in this particular representation, which ignores the Hund 
rules and, therefore, ignores an essential part of the problem. Moreover, even in this representation the degeneracy 
must enter the equations explicitly, since it always exist if we ignore the intra shell spin-spin exchange (first Hund 
rule) and orbital-orbital exchange (second Hund rule). For the less than half filled shell this is not very essential 
unless we consider spin excitations and do not take into account the crystal field since the requirement of maximal 
spin and orbital moment makes the ground-state wave function unique. Let us rewrite our equations in the following 
matrix form, taking into account the above assumption (here, for a brevity we omit the k-dependence in /i(k), V^(k), 
and 0(k)): 



^ Without coupling to the lattice this state may be degenerate. 
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The on-site matrix elements {f^)"'"^, which in the case of an orthogonal basis set give "selection rules", give a 
non-zero value: {flY"'^ — (S^\fwi • fl\^) = '^Ziy For example, for the case of one c-orbital and two delocalized /- 
orbitals Q-^^ = [(1 - 5uy^)0^cOcv, +5yvAOcc - |0,.cP)]/ det O, where detO = Occ - |C^cP - Py.A'^- If we restrict 
ourselves with the accuracy ^ o{0^), i.e. neglect the matrix elements ~ lOi/cP, then 0'^^^, ~ S^^u. If we ignore the 
multiplet structure, then all energies = A2 coincide, and we find that the combination of population numbers 
p[r„,r„+i] _ S'r„+i,r„+i^r„ + 5'r„,r„^r„+i for different transitions, a, does not depend on a: P"^ — P ^. For the 
states r„ we have '5r„.r„ = (r„|r„) = (0|0) = 1 and for r„+i 5r„+ij„+i = (r„+i|r„+i> - (6|/,/t|0) ~ 1. With 
these approximations we obtain the following system of equations 

[&LL.Ju;-Ol\hl^^^-Ollv:^P] Pi]A2+OiIVl,A \(Gfl,{Kzu) Gif:,(k,za;)\ 



OA' oilp 



(20) 



0-l,P 5,.,P 

Let us now introduce the effective Fermions / = X/ \fP with the help of the following transformation: 

(Ia;-fl)G = P ^ U(Iw - n)U"^g = UPU, (21) 

where the block-matrix U is: 

U.(;^) (22, 



and g ~ UGU. After this transformation Eq.(20) acquires the following form 



{&LL.iu:~oi\hi^^^-{oil^){v:^^^)\ [(oz_i,VP)A2 + ozi,(^L,.VP)] \ / atA^^^^) Qlt'i^M 



Let us compare the system of equations, Eq.(|23[), for the GFs, G, written in orbital representation, with the system 
of equations for the Fermion GFs, F, corresponding to the secular problem in normal electronic structure calculations, 
Eq.([7|). We find that the systems coincide if we renormalize mixing and overlap matrices as follows: 

OA ^ OA = OA Vl., - Vl., - Vl., Vp, (24) 

and put 0~^, — 5^^'- Thus, the matrix equation for the renormalized GFs, (Iw — O^^Y)^ = O^^, can be rewritten 
in the form containing the standard secular matrix used in the electronic structure calculations 

(0^-V)g = L (25) 

Now, using the standard procedure of Cholesky decomposition, I ={Ou) — V)g = (ZwZ — V)g, we can rewrite this 
equation in terms of an effective Hamiltonian: 

(Iw - Z"^VZ"^)ZgZ = I, (26) 

or, 

(I^ - h)g = L (27) 

Thus, diagonalizing this equation in k-space, we find the renormalized band structure. The parameter P, is calculated 
from the Green function G^^ . 

In the paper IlQ we have shown that within the ARF each graph contributing to the Green function and the self- 
energy in the weak coupling perturbation theory has its counter-part in the strong-coupling theory and, therefore, 
one can use the same expression for the total energy in the ARF as in the case of weak-coupling theory but with the 
renormalizations shown in Eq.(p3). 



^ In general, of course, this is a matrix M x M, where M is number of transitions. 
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III. ELECTRONIC STRUCTURE CALCULATION IN THE ARE 



Since the secular equations to be solved include, in an explicit way, many body corrections in ARF, we will refer 
to this calculation as an electronic structure calculation in ARF (ESC- ARF). Let us now describe the numerical 
implementation of the ESC-ARF method. Following the analysis of the previous section, as well as our previous 
works, the f-states are divided into a lower level and an upper level. In ESC the electrons are normally separated 
into core and valence electrons, and it is not hard, in practice, to include the many body corrections derived in the 
previous section. The lower transitions are simply single electron core levels in the ESC, and the upper transitions 
are all described by the set of equations described in the previous section. Due to the analogy made between the 
obtained equations, Eq.([23|) and the standard secular equation solved for the Kohn-Sham equation, Eq.(^, one needs 
to perform only relatively small modifications to any electronic structure method. Implemented in a full potential, 
site centered basis function method (to be specific, a linear muffin-tin orbital method) the Hamiltonian matrix has 
the schematic form as is shown in Fig. |l|. The energy of the upper transition, A2 may be used as a variable input 
parameter, but it can also be taken from bremsstrahlung isochromat spectra (BIS). The values obtained in this way 
are presented in Table |. 



element 


Co 


Pr 


Nd 




Sm 


/-peak 


3.52 


2.27 


1.68 


1.08'' 


0.64 



TABLE I. Positions of the /-peak above the Fermi level in the lanthanide series in eV, taken from BIS, from ren 
^Pm is radioactive, this value is a fit between Nd and Sm 



The renormalizing factor, v P, which should multiply the Hamiltonian and the overlap matrices in the places 
specified in Fig. |l| has to be calculated sclf-consistently. Since we only take into account the states, r„ and r„+i, 
there are only two non-zero population numbers, 7Vr„ and A^r„+u which gives rise to a single P. Further, we disregard 
the multiplet structure and therefore all transitions have the same energy, Ai for the lower level and A2 from the 
upper level. The population numbers distribute as follows; 

Hf = n- Ni + {n + l)J2,,N2'^^^ n + T] /2o^ 

where n and n + I are the degeneracy of the states r„ and r„+i, respectively, and v is an index of the available 
single-electron orbitals to construct r„_|_i from r„. If we assume that all upper transitions (this number is the same 
as the number of available single-electron orbitals) are occupied symmetrically, we find that X]jy^2 = "-^21 i-G- 
-^2 = ^7/(14 — n) and A^i = 1 — 77. Then, the renormalizing factor will be 



14- 



The /-electron locator GF becomes 



and 



lUJn - Ai lUn - A2 

-, _ 14-(ri.-|-l) 14-(»+l) 

14-n ^ , 14-Ti ^ /oi\ 

y<^ = — ■■ T ^ 7—- (-^1) 

lUJn - A2 lUJn - A3 

Here, A3 is the energy of the n+1 to n + 2 transition, which is much above the Fermi level and can be ignored. As 
seen the equations above contain an unknown factor rj. This factor should be calculated from the Green function 

Tl = --y Im{g,{Lo + i5,'k)f{Lo)). (32) 

This equation has to be solved self-consistently, since the factor P renormalizes the Hamiltonian and^pverlap. 

The corrections has been implemented in a full potential linear muffin tin (FP-LMTO) methodQa. For the lan- 
thanides, the face centered cubic structure (fee) cubic lattice was used. This was due to the simplicity to perform 
calculations with one atom per primitive cell. In our calculations we set the same criteria for convergence in k-point 
summation, self-consistency, basis set truncation etc. as is done in normal electronic structure calculations. 
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IV. RESULTS 



We selected Ce, Pr, Nd and Sm as elements of special interest to test our theory for, since we, in a previous work 
showed that the many-body effect of the upper transition is largest for them. The renormalizing factor, a/P, is, 
according to our calculations, less than unity, and, thus, the mixing matrix elements are reduced somewhat. This 
means that the bandwidth is reduced compared to that obtained from a regular ESC. The upper transition can be 
taken from experiment, when available, or alternatively calculated from a constrained total energy calculation. In 
the present investigation we also made a scan, changing A2, from a position high above the Fermi- level, to a place 
closer to the Fermi-level. The results of the ESC-ARF theory are shown in Fig. g. In this figure, the top figure shows 
the volumes as a function of the position of the upper transition, A2. The central one, shows the f-occupation and 
the lower one shows the renormalization factor, \/P, for Ce, Pr, Nd and Pm. We note that the renormalizing factor 
is quite close to unity for all cases which means that the effect of the correlations for these materials is quite small. 
Further, when the upper transition is lowered, the /-occupation increases, and there is a contribution to the cohesive 
energy from the /-electrons, that results in that the volume decreases. In Table || we list the volumes for the selected 
lanthanides. 



element Expt.'" LDA*" here 

Co 3.8115 3.8994 3.7485 

Pr 3.8178 3.8375 3.9030 

Nd 3.8045 3.7848 3.7484 

Pm'^ 3.7832 3.7415 3.5182 



TABLE IL Volumes, Rsws, for the light lanthanides, in a.u. When the position of the upper transition, A2 is put to the 
experimental value shown in Table | 
''Taken from refE. 

''ESC with the 4/ treated as core electrons. 

"^Pm is radioactive, this value is a fit between Nd and Sm. 

In Fig. ^, the density of states after the self-consistent cycle is shown. The occupied part of the d-band is quite 
similar for the case of a normal calculation treating the f-level as part of the core. Also, the 4/-band like feature 
is rather similar to what would be obtained from a regular ESC calculation, with a small narrowing due to the 
renormalization and a shift to higher energies, so that the number of f-electrons occupying this feature is small, ~ 
0.05-0.1. The resulting electronic structure of our ESC-ARF theory is actually rather similar to that obtained when 
including Hubbard-Kanamory interactions to an LDA Hamiltonian (or the related LDA+U method) and to what is 
obtained from calculations including self-interaction corrections. 

In Fig. I the f-bandwidth as a function of the position of the upper transition with respect to the Fermi level is 
plotted. The circles correspond to a self-consistent calculation where yfP is put to 1 during the whole cycle. The 
boxes are the result from a full calculation including all renormalization effects described above. A band narrowing 
effect is seen when \fP 7^ 1, but one may note that the effect is not extremely pronounced for the presently studied 
elements. One may speculate that for other systems, where a larger population of the upper level is expected (say, 
the level A2 is much closer to Fermi enp^y), the renormalization is more pronounced. Since the present modifications 
are similar to those presented for (5-Pul3, that does show a larger population of the upper band, it may very well be 
so that the present corrections are most important for this element. 



V. CONCLUSION 



We have presented a method to incorporate strong correlations into electronic structure calculations (ESC-ARF) 
and performed numerical tests of the theory. The results show a narrowing of the 4/-bands compared to LDA (or 
GGA) results. Due to the division of the 4/-states in a higher and lower level, we find that the cohesion from the 
4/-states is non-zero (in contrast to calculations that treat the 4/-states as core states) but substantially smaller than 
what it would have been in an LDA (or GGA) calculation. 

The present results are rather similar to our previous calculations, that include modifications in an LMTO-ASA 
method, although the full-potential treatment does include finer details of the potential. Due to the generalization of 
the corrections to an electronic structure method, that can treat any crystal structure, we believe the present method 
should have a potential in calculating the electronic structure of /-based materials that are hard to treat in a regular 
LDA (or GGA) ESC methods. 
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FIG. 1. The Hamiltonian (and overlap) matrix modified according to the text. A2 is the position of the upper transition. 
FIG. 2. The volumes, f-occupation and \/P-factor for Ce, Pr, Nd and Pm, as a function of the position of the upper transition, 



A. 



2- 



FIG. 3. The d- and /-projected density of states for the two cases when there are f-states in the valence basis with all 
corrections described above, and when a normal ESC with the /-electrons treated as core. 

FIG. 4. The bandwidth of the /-states, when the renormalizing factor is put to one and when it is not. 
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